\(\int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx\) [265]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (verified)
   Maple [B] (warning: unable to verify)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [B] (verification not implemented)
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 35, antiderivative size = 156 \[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\frac {(19 A+5 B) \text {arctanh}\left (\frac {\sqrt {a} \sqrt {\sec (c+d x)} \sin (c+d x)}{\sqrt {2} \sqrt {a+a \sec (c+d x)}}\right )}{16 \sqrt {2} a^{5/2} d}-\frac {(A-B) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{4 d (a+a \sec (c+d x))^{5/2}}-\frac {(9 A-B) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{16 a d (a+a \sec (c+d x))^{3/2}} \]

[Out]

-1/4*(A-B)*sec(d*x+c)^(3/2)*sin(d*x+c)/d/(a+a*sec(d*x+c))^(5/2)-1/16*(9*A-B)*sec(d*x+c)^(3/2)*sin(d*x+c)/a/d/(
a+a*sec(d*x+c))^(3/2)+1/32*(19*A+5*B)*arctanh(1/2*sin(d*x+c)*a^(1/2)*sec(d*x+c)^(1/2)*2^(1/2)/(a+a*sec(d*x+c))
^(1/2))/a^(5/2)/d*2^(1/2)

Rubi [A] (verified)

Time = 0.72 (sec) , antiderivative size = 203, normalized size of antiderivative = 1.30, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {4104, 4105, 4098, 3893, 212} \[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\frac {(19 A+5 B) \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x) \sqrt {\sec (c+d x)}}{\sqrt {2} \sqrt {a \sec (c+d x)+a}}\right )}{16 \sqrt {2} a^{5/2} d}-\frac {(9 A-B) \sin (c+d x) \sqrt {\sec (c+d x)}}{16 a^2 d \sqrt {a \sec (c+d x)+a}}+\frac {(5 A+3 B) \sin (c+d x) \sqrt {\sec (c+d x)}}{16 a d (a \sec (c+d x)+a)^{3/2}}+\frac {(A-B) \sin (c+d x) \sqrt {\sec (c+d x)}}{4 d (a \sec (c+d x)+a)^{5/2}} \]

[In]

Int[(Sqrt[Sec[c + d*x]]*(A + B*Sec[c + d*x]))/(a + a*Sec[c + d*x])^(5/2),x]

[Out]

((19*A + 5*B)*ArcTanh[(Sqrt[a]*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(Sqrt[2]*Sqrt[a + a*Sec[c + d*x]])])/(16*Sqrt[
2]*a^(5/2)*d) + ((A - B)*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(4*d*(a + a*Sec[c + d*x])^(5/2)) + ((5*A + 3*B)*Sqrt
[Sec[c + d*x]]*Sin[c + d*x])/(16*a*d*(a + a*Sec[c + d*x])^(3/2)) - ((9*A - B)*Sqrt[Sec[c + d*x]]*Sin[c + d*x])
/(16*a^2*d*Sqrt[a + a*Sec[c + d*x]])

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 3893

Int[Sqrt[csc[(e_.) + (f_.)*(x_)]*(d_.)]/Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Dist[-2*b*(d/
(a*f)), Subst[Int[1/(2*b - d*x^2), x], x, b*(Cot[e + f*x]/(Sqrt[a + b*Csc[e + f*x]]*Sqrt[d*Csc[e + f*x]]))], x
] /; FreeQ[{a, b, d, e, f}, x] && EqQ[a^2 - b^2, 0]

Rule 4098

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> Simp[A*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*((d*Csc[e + f*x])^n/(f*n)), x] - Dist[(
a*A*m - b*B*n)/(b*d*n), Int[(a + b*Csc[e + f*x])^m*(d*Csc[e + f*x])^(n + 1), x], x] /; FreeQ[{a, b, d, e, f, A
, B, m, n}, x] && NeQ[A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && EqQ[m + n + 1, 0] &&  !LeQ[m, -1]

Rule 4104

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> Simp[d*(A*b - a*B)*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*((d*Csc[e + f*x])^(n - 1)/(
a*f*(2*m + 1))), x] - Dist[1/(a*b*(2*m + 1)), Int[(a + b*Csc[e + f*x])^(m + 1)*(d*Csc[e + f*x])^(n - 1)*Simp[A
*(a*d*(n - 1)) - B*(b*d*(n - 1)) - d*(a*B*(m - n + 1) + A*b*(m + n))*Csc[e + f*x], x], x], x] /; FreeQ[{a, b,
d, e, f, A, B}, x] && NeQ[A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && LtQ[m, -2^(-1)] && GtQ[n, 0]

Rule 4105

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> Simp[(-(A*b - a*B))*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*((d*Csc[e + f*x])^n/(b*f*(
2*m + 1))), x] - Dist[1/(a^2*(2*m + 1)), Int[(a + b*Csc[e + f*x])^(m + 1)*(d*Csc[e + f*x])^n*Simp[b*B*n - a*A*
(2*m + n + 1) + (A*b - a*B)*(m + n + 1)*Csc[e + f*x], x], x], x] /; FreeQ[{a, b, d, e, f, A, B, n}, x] && NeQ[
A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && LtQ[m, -2^(-1)] &&  !GtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {(A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{4 d (a+a \sec (c+d x))^{5/2}}+\frac {\int \frac {-\frac {1}{2} a (A-B)+2 a (A+B) \sec (c+d x)}{\sqrt {\sec (c+d x)} (a+a \sec (c+d x))^{3/2}} \, dx}{4 a^2} \\ & = \frac {(A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{4 d (a+a \sec (c+d x))^{5/2}}+\frac {(5 A+3 B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a d (a+a \sec (c+d x))^{3/2}}+\frac {\int \frac {-\frac {1}{4} a^2 (9 A-B)+\frac {1}{2} a^2 (5 A+3 B) \sec (c+d x)}{\sqrt {\sec (c+d x)} \sqrt {a+a \sec (c+d x)}} \, dx}{8 a^4} \\ & = \frac {(A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{4 d (a+a \sec (c+d x))^{5/2}}+\frac {(5 A+3 B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a d (a+a \sec (c+d x))^{3/2}}-\frac {(9 A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a^2 d \sqrt {a+a \sec (c+d x)}}+\frac {(19 A+5 B) \int \frac {\sqrt {\sec (c+d x)}}{\sqrt {a+a \sec (c+d x)}} \, dx}{32 a^2} \\ & = \frac {(A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{4 d (a+a \sec (c+d x))^{5/2}}+\frac {(5 A+3 B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a d (a+a \sec (c+d x))^{3/2}}-\frac {(9 A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a^2 d \sqrt {a+a \sec (c+d x)}}-\frac {(19 A+5 B) \text {Subst}\left (\int \frac {1}{2 a-x^2} \, dx,x,-\frac {a \sqrt {\sec (c+d x)} \sin (c+d x)}{\sqrt {a+a \sec (c+d x)}}\right )}{16 a^2 d} \\ & = \frac {(19 A+5 B) \text {arctanh}\left (\frac {\sqrt {a} \sqrt {\sec (c+d x)} \sin (c+d x)}{\sqrt {2} \sqrt {a+a \sec (c+d x)}}\right )}{16 \sqrt {2} a^{5/2} d}+\frac {(A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{4 d (a+a \sec (c+d x))^{5/2}}+\frac {(5 A+3 B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a d (a+a \sec (c+d x))^{3/2}}-\frac {(9 A-B) \sqrt {\sec (c+d x)} \sin (c+d x)}{16 a^2 d \sqrt {a+a \sec (c+d x)}} \\ \end{align*}

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(443\) vs. \(2(156)=312\).

Time = 6.10 (sec) , antiderivative size = 443, normalized size of antiderivative = 2.84 \[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=-\frac {A \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{4 d (a (1+\sec (c+d x)))^{5/2}}-\frac {B \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{4 d (a (1+\sec (c+d x)))^{5/2}}-\frac {9 A \sec ^{\frac {3}{2}}(c+d x) (1+\sec (c+d x)) \sin (c+d x)}{16 d (a (1+\sec (c+d x)))^{5/2}}-\frac {19 A \arctan \left (\frac {\sqrt {2} \sqrt {\sec (c+d x)}}{\sqrt {1-\sec (c+d x)}}\right ) (1+\sec (c+d x))^2 \tan (c+d x)}{16 \sqrt {2} d \sqrt {1-\sec (c+d x)} (a (1+\sec (c+d x)))^{5/2}}-\frac {5 B (1+\sec (c+d x))^{5/2} \left (\frac {2 \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{d (1+\sec (c+d x))^{3/2}}-\frac {\left (2 \arcsin \left (\sqrt {\sec (c+d x)}\right )-\sqrt {2} \arctan \left (\frac {\sqrt {2} \sqrt {\sec (c+d x)}}{\sqrt {1-\sec (c+d x)}}\right )\right ) \tan (c+d x)}{d \sqrt {1-\sec (c+d x)} \sqrt {1+\sec (c+d x)}}-\frac {2 \left (\arcsin \left (\sqrt {1-\sec (c+d x)}\right )+\sqrt {1-\sec (c+d x)} \sqrt {\sec (c+d x)}\right ) \tan (c+d x)}{d \sqrt {1-\sec (c+d x)} \sqrt {1+\sec (c+d x)}}\right )}{32 (a (1+\sec (c+d x)))^{5/2}} \]

[In]

Integrate[(Sqrt[Sec[c + d*x]]*(A + B*Sec[c + d*x]))/(a + a*Sec[c + d*x])^(5/2),x]

[Out]

-1/4*(A*Sec[c + d*x]^(3/2)*Sin[c + d*x])/(d*(a*(1 + Sec[c + d*x]))^(5/2)) - (B*Sec[c + d*x]^(5/2)*Sin[c + d*x]
)/(4*d*(a*(1 + Sec[c + d*x]))^(5/2)) - (9*A*Sec[c + d*x]^(3/2)*(1 + Sec[c + d*x])*Sin[c + d*x])/(16*d*(a*(1 +
Sec[c + d*x]))^(5/2)) - (19*A*ArcTan[(Sqrt[2]*Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]]*(1 + Sec[c + d*x])^2
*Tan[c + d*x])/(16*Sqrt[2]*d*Sqrt[1 - Sec[c + d*x]]*(a*(1 + Sec[c + d*x]))^(5/2)) - (5*B*(1 + Sec[c + d*x])^(5
/2)*((2*Sec[c + d*x]^(5/2)*Sin[c + d*x])/(d*(1 + Sec[c + d*x])^(3/2)) - ((2*ArcSin[Sqrt[Sec[c + d*x]]] - Sqrt[
2]*ArcTan[(Sqrt[2]*Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]])*Tan[c + d*x])/(d*Sqrt[1 - Sec[c + d*x]]*Sqrt[1
 + Sec[c + d*x]]) - (2*(ArcSin[Sqrt[1 - Sec[c + d*x]]] + Sqrt[1 - Sec[c + d*x]]*Sqrt[Sec[c + d*x]])*Tan[c + d*
x])/(d*Sqrt[1 - Sec[c + d*x]]*Sqrt[1 + Sec[c + d*x]])))/(32*(a*(1 + Sec[c + d*x]))^(5/2))

Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(418\) vs. \(2(131)=262\).

Time = 5.01 (sec) , antiderivative size = 419, normalized size of antiderivative = 2.69

method result size
default \(-\frac {\sqrt {-\frac {\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}+1}{\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\, \left (\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1\right ) \sqrt {-\frac {2 a}{\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\, \left (2 A \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (1-\cos \left (d x +c \right )\right )^{3} \csc \left (d x +c \right )^{3}-2 B \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (1-\cos \left (d x +c \right )\right )^{3} \csc \left (d x +c \right )^{3}-11 A \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )+3 B \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )-19 A \arctan \left (\frac {-\cot \left (d x +c \right )+\csc \left (d x +c \right )}{\sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\right )-5 B \arctan \left (\frac {-\cot \left (d x +c \right )+\csc \left (d x +c \right )}{\sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\right )\right )}{32 a^{3} d \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\) \(419\)
parts \(-\frac {A \sqrt {-\frac {\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}+1}{\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\, \left (\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1\right ) \sqrt {-\frac {2 a}{\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\, \left (2 \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (1-\cos \left (d x +c \right )\right )^{3} \csc \left (d x +c \right )^{3}-11 \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )-19 \arctan \left (\frac {-\cot \left (d x +c \right )+\csc \left (d x +c \right )}{\sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\right )\right )}{32 d \,a^{3} \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}-\frac {B \left (-\frac {\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}+1}{\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\right )^{\frac {3}{2}} \left (\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1\right )^{2} \sqrt {-\frac {2 a}{\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\, \left (2 \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (1-\cos \left (d x +c \right )\right )^{3} \csc \left (d x +c \right )^{3}-3 \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}\, \left (-\cot \left (d x +c \right )+\csc \left (d x +c \right )\right )+5 \arctan \left (\frac {-\cot \left (d x +c \right )+\csc \left (d x +c \right )}{\sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\right )\right )}{32 d \,a^{3} \left (\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}+1\right ) \sqrt {-\left (1-\cos \left (d x +c \right )\right )^{2} \csc \left (d x +c \right )^{2}-1}}\) \(583\)

[In]

int((A+B*sec(d*x+c))*sec(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(5/2),x,method=_RETURNVERBOSE)

[Out]

-1/32/a^3/d*(-((1-cos(d*x+c))^2*csc(d*x+c)^2+1)/((1-cos(d*x+c))^2*csc(d*x+c)^2-1))^(1/2)*((1-cos(d*x+c))^2*csc
(d*x+c)^2-1)*(-2*a/((1-cos(d*x+c))^2*csc(d*x+c)^2-1))^(1/2)*(2*A*(-(1-cos(d*x+c))^2*csc(d*x+c)^2-1)^(1/2)*(1-c
os(d*x+c))^3*csc(d*x+c)^3-2*B*(-(1-cos(d*x+c))^2*csc(d*x+c)^2-1)^(1/2)*(1-cos(d*x+c))^3*csc(d*x+c)^3-11*A*(-(1
-cos(d*x+c))^2*csc(d*x+c)^2-1)^(1/2)*(-cot(d*x+c)+csc(d*x+c))+3*B*(-(1-cos(d*x+c))^2*csc(d*x+c)^2-1)^(1/2)*(-c
ot(d*x+c)+csc(d*x+c))-19*A*arctan(1/(-(1-cos(d*x+c))^2*csc(d*x+c)^2-1)^(1/2)*(-cot(d*x+c)+csc(d*x+c)))-5*B*arc
tan(1/(-(1-cos(d*x+c))^2*csc(d*x+c)^2-1)^(1/2)*(-cot(d*x+c)+csc(d*x+c))))/(-(1-cos(d*x+c))^2*csc(d*x+c)^2-1)^(
1/2)

Fricas [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 502, normalized size of antiderivative = 3.22 \[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\left [\frac {\sqrt {2} {\left ({\left (19 \, A + 5 \, B\right )} \cos \left (d x + c\right )^{3} + 3 \, {\left (19 \, A + 5 \, B\right )} \cos \left (d x + c\right )^{2} + 3 \, {\left (19 \, A + 5 \, B\right )} \cos \left (d x + c\right ) + 19 \, A + 5 \, B\right )} \sqrt {a} \log \left (-\frac {a \cos \left (d x + c\right )^{2} - 2 \, \sqrt {2} \sqrt {a} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right ) - 2 \, a \cos \left (d x + c\right ) - 3 \, a}{\cos \left (d x + c\right )^{2} + 2 \, \cos \left (d x + c\right ) + 1}\right ) - \frac {4 \, {\left ({\left (13 \, A - 5 \, B\right )} \cos \left (d x + c\right )^{2} + {\left (9 \, A - B\right )} \cos \left (d x + c\right )\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sin \left (d x + c\right )}{\sqrt {\cos \left (d x + c\right )}}}{64 \, {\left (a^{3} d \cos \left (d x + c\right )^{3} + 3 \, a^{3} d \cos \left (d x + c\right )^{2} + 3 \, a^{3} d \cos \left (d x + c\right ) + a^{3} d\right )}}, -\frac {\sqrt {2} {\left ({\left (19 \, A + 5 \, B\right )} \cos \left (d x + c\right )^{3} + 3 \, {\left (19 \, A + 5 \, B\right )} \cos \left (d x + c\right )^{2} + 3 \, {\left (19 \, A + 5 \, B\right )} \cos \left (d x + c\right ) + 19 \, A + 5 \, B\right )} \sqrt {-a} \arctan \left (\frac {\sqrt {2} \sqrt {-a} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt {\cos \left (d x + c\right )}}{a \sin \left (d x + c\right )}\right ) + \frac {2 \, {\left ({\left (13 \, A - 5 \, B\right )} \cos \left (d x + c\right )^{2} + {\left (9 \, A - B\right )} \cos \left (d x + c\right )\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sin \left (d x + c\right )}{\sqrt {\cos \left (d x + c\right )}}}{32 \, {\left (a^{3} d \cos \left (d x + c\right )^{3} + 3 \, a^{3} d \cos \left (d x + c\right )^{2} + 3 \, a^{3} d \cos \left (d x + c\right ) + a^{3} d\right )}}\right ] \]

[In]

integrate((A+B*sec(d*x+c))*sec(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(5/2),x, algorithm="fricas")

[Out]

[1/64*(sqrt(2)*((19*A + 5*B)*cos(d*x + c)^3 + 3*(19*A + 5*B)*cos(d*x + c)^2 + 3*(19*A + 5*B)*cos(d*x + c) + 19
*A + 5*B)*sqrt(a)*log(-(a*cos(d*x + c)^2 - 2*sqrt(2)*sqrt(a)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sqrt(cos(
d*x + c))*sin(d*x + c) - 2*a*cos(d*x + c) - 3*a)/(cos(d*x + c)^2 + 2*cos(d*x + c) + 1)) - 4*((13*A - 5*B)*cos(
d*x + c)^2 + (9*A - B)*cos(d*x + c))*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sin(d*x + c)/sqrt(cos(d*x + c)))/
(a^3*d*cos(d*x + c)^3 + 3*a^3*d*cos(d*x + c)^2 + 3*a^3*d*cos(d*x + c) + a^3*d), -1/32*(sqrt(2)*((19*A + 5*B)*c
os(d*x + c)^3 + 3*(19*A + 5*B)*cos(d*x + c)^2 + 3*(19*A + 5*B)*cos(d*x + c) + 19*A + 5*B)*sqrt(-a)*arctan(sqrt
(2)*sqrt(-a)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))/(a*sin(d*x + c))) + 2*((13*A - 5*B)*co
s(d*x + c)^2 + (9*A - B)*cos(d*x + c))*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sin(d*x + c)/sqrt(cos(d*x + c))
)/(a^3*d*cos(d*x + c)^3 + 3*a^3*d*cos(d*x + c)^2 + 3*a^3*d*cos(d*x + c) + a^3*d)]

Sympy [F]

\[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\int \frac {\left (A + B \sec {\left (c + d x \right )}\right ) \sqrt {\sec {\left (c + d x \right )}}}{\left (a \left (\sec {\left (c + d x \right )} + 1\right )\right )^{\frac {5}{2}}}\, dx \]

[In]

integrate((A+B*sec(d*x+c))*sec(d*x+c)**(1/2)/(a+a*sec(d*x+c))**(5/2),x)

[Out]

Integral((A + B*sec(c + d*x))*sqrt(sec(c + d*x))/(a*(sec(c + d*x) + 1))**(5/2), x)

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 5924 vs. \(2 (131) = 262\).

Time = 0.91 (sec) , antiderivative size = 5924, normalized size of antiderivative = 37.97 \[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\text {Too large to display} \]

[In]

integrate((A+B*sec(d*x+c))*sec(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(5/2),x, algorithm="maxima")

[Out]

1/32*((19*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x
 + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(4*d*x + 4*c)^2 + 304*(log(cos(1/2*d*x
+ 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x +
 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(3*d*x + 3*c)^2 + 684*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x +
1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x +
 1/2*c) + 1))*cos(2*d*x + 2*c)^2 + 304*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x +
1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(d*x + c)^
2 + 19*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x +
1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(4*d*x + 4*c)^2 + 304*(log(cos(1/2*d*x + 1
/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/
2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(3*d*x + 3*c)^2 + 684*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2
*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/
2*c) + 1))*sin(2*d*x + 2*c)^2 + 304*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2
*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(d*x + c)^2 +
 2*(76*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x +
1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(3*d*x + 3*c) + 114*(log(cos(1/2*d*x + 1/2
*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*
c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(2*d*x + 2*c) + 76*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2
 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c)
+ 1))*cos(d*x + c) + 19*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 19
*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 26*sin(7/2*d*x + 7/2*c) -
 10*sin(5/2*d*x + 5/2*c) + 10*sin(3/2*d*x + 3/2*c) + 26*sin(1/2*d*x + 1/2*c))*cos(4*d*x + 4*c) + 104*(2*sin(3*
d*x + 3*c) + 3*sin(2*d*x + 2*c) + 2*sin(d*x + c))*cos(7/2*d*x + 7/2*c) + 8*(114*(log(cos(1/2*d*x + 1/2*c)^2 +
sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2
*sin(1/2*d*x + 1/2*c) + 1))*cos(2*d*x + 2*c) + 76*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin
(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*co
s(d*x + c) + 19*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 19*log(cos
(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 10*sin(5/2*d*x + 5/2*c) + 10*sin(
3/2*d*x + 3/2*c) + 26*sin(1/2*d*x + 1/2*c))*cos(3*d*x + 3*c) + 40*(3*sin(2*d*x + 2*c) + 2*sin(d*x + c))*cos(5/
2*d*x + 5/2*c) + 12*(76*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - l
og(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(d*x + c) + 19*log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 19*log(cos(1/2*d*x + 1/2*c)^2 + sin(
1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) + 10*sin(3/2*d*x + 3/2*c) + 26*sin(1/2*d*x + 1/2*c))*cos(2*d*
x + 2*c) + 8*(19*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 19*log(co
s(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) + 26*sin(1/2*d*x + 1/2*c))*cos(d*x
 + c) + 4*(38*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(3*d*x + 3*c) + 57*(log(cos(1/2*d*x
 + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x
+ 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(2*d*x + 2*c) + 38*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/
2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1
/2*c) + 1))*sin(d*x + c) + 13*cos(7/2*d*x + 7/2*c) + 5*cos(5/2*d*x + 5/2*c) - 5*cos(3/2*d*x + 3/2*c) - 13*cos(
1/2*d*x + 1/2*c))*sin(4*d*x + 4*c) - 52*(4*cos(3*d*x + 3*c) + 6*cos(2*d*x + 2*c) + 4*cos(d*x + c) + 1)*sin(7/2
*d*x + 7/2*c) + 16*(57*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - lo
g(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(2*d*x + 2*c) + 38*(log(co
s(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin
(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(d*x + c) + 5*cos(5/2*d*x + 5/2*c) - 5*cos(3/2*d*x + 3/2
*c) - 13*cos(1/2*d*x + 1/2*c))*sin(3*d*x + 3*c) - 20*(6*cos(2*d*x + 2*c) + 4*cos(d*x + c) + 1)*sin(5/2*d*x + 5
/2*c) + 24*(38*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(d*x + c) - 5*cos(3/2*d*x + 3/2*c)
 - 13*cos(1/2*d*x + 1/2*c))*sin(2*d*x + 2*c) + 20*(4*cos(d*x + c) + 1)*sin(3/2*d*x + 3/2*c) - 80*cos(3/2*d*x +
 3/2*c)*sin(d*x + c) - 208*cos(1/2*d*x + 1/2*c)*sin(d*x + c) + 19*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1
/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 19*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x
 + 1/2*c) + 1) + 52*sin(1/2*d*x + 1/2*c))*A/((sqrt(2)*a^2*cos(4*d*x + 4*c)^2 + 16*sqrt(2)*a^2*cos(3*d*x + 3*c)
^2 + 36*sqrt(2)*a^2*cos(2*d*x + 2*c)^2 + 16*sqrt(2)*a^2*cos(d*x + c)^2 + sqrt(2)*a^2*sin(4*d*x + 4*c)^2 + 16*s
qrt(2)*a^2*sin(3*d*x + 3*c)^2 + 36*sqrt(2)*a^2*sin(2*d*x + 2*c)^2 + 48*sqrt(2)*a^2*sin(2*d*x + 2*c)*sin(d*x +
c) + 16*sqrt(2)*a^2*sin(d*x + c)^2 + 8*sqrt(2)*a^2*cos(d*x + c) + sqrt(2)*a^2 + 2*(4*sqrt(2)*a^2*cos(3*d*x + 3
*c) + 6*sqrt(2)*a^2*cos(2*d*x + 2*c) + 4*sqrt(2)*a^2*cos(d*x + c) + sqrt(2)*a^2)*cos(4*d*x + 4*c) + 8*(6*sqrt(
2)*a^2*cos(2*d*x + 2*c) + 4*sqrt(2)*a^2*cos(d*x + c) + sqrt(2)*a^2)*cos(3*d*x + 3*c) + 12*(4*sqrt(2)*a^2*cos(d
*x + c) + sqrt(2)*a^2)*cos(2*d*x + 2*c) + 4*(2*sqrt(2)*a^2*sin(3*d*x + 3*c) + 3*sqrt(2)*a^2*sin(2*d*x + 2*c) +
 2*sqrt(2)*a^2*sin(d*x + c))*sin(4*d*x + 4*c) + 16*(3*sqrt(2)*a^2*sin(2*d*x + 2*c) + 2*sqrt(2)*a^2*sin(d*x + c
))*sin(3*d*x + 3*c))*sqrt(a)) + (4*(3*sin(3/2*d*x + 3/2*c) + 5*sin(7/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d
*x + 3/2*c))) - 3*sin(5/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 5*sin(1/3*arctan2(sin(3/2*d*x
 + 3/2*c), cos(3/2*d*x + 3/2*c))))*cos(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 40*(2*sin(3*
d*x + 3*c) + 3*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 2*sin(2/3*arctan2(sin(3/2*d*x +
3/2*c), cos(3/2*d*x + 3/2*c))))*cos(7/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 24*(2*sin(3*d*x
 + 3*c) + 3*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 2*sin(2/3*arctan2(sin(3/2*d*x + 3/2
*c), cos(3/2*d*x + 3/2*c))))*cos(5/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 24*(3*sin(3/2*d*x
+ 3/2*c) - 5*sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c
), cos(3/2*d*x + 3/2*c))) + 16*(3*sin(3/2*d*x + 3/2*c) - 5*sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x +
 3/2*c))))*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 5*(16*cos(3*d*x + 3*c)^2 + 2*(4*cos(
3*d*x + 3*c) + 6*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 4*cos(2/3*arctan2(sin(3/2*d*x
+ 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*cos(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + cos(8/3*a
rctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 12*(4*cos(3*d*x + 3*c) + 4*cos(2/3*arctan2(sin(3/2*d*x
 + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 36*cos(4
/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 8*(4*cos(3*d*x + 3*c) + 1)*cos(2/3*arctan2(sin(3/2
*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 16*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 16
*sin(3*d*x + 3*c)^2 + 4*(2*sin(3*d*x + 3*c) + 3*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) +
 2*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2
*d*x + 3/2*c))) + sin(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 48*(sin(3*d*x + 3*c) + sin(
2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x +
3/2*c))) + 36*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 32*sin(3*d*x + 3*c)*sin(2/3*arc
tan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 16*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2
*c)))^2 + 8*cos(3*d*x + 3*c) + 1)*log(cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + sin(1/3
*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 2*sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x
+ 3/2*c))) + 1) - 5*(16*cos(3*d*x + 3*c)^2 + 2*(4*cos(3*d*x + 3*c) + 6*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), c
os(3/2*d*x + 3/2*c))) + 4*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*cos(8/3*arctan2(si
n(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + cos(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 +
12*(4*cos(3*d*x + 3*c) + 4*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*cos(4/3*arctan2(s
in(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 36*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^
2 + 8*(4*cos(3*d*x + 3*c) + 1)*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 16*cos(2/3*arcta
n2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 16*sin(3*d*x + 3*c)^2 + 4*(2*sin(3*d*x + 3*c) + 3*sin(4/3*
arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 2*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3
/2*c))))*sin(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + sin(8/3*arctan2(sin(3/2*d*x + 3/2*c),
cos(3/2*d*x + 3/2*c)))^2 + 48*(sin(3*d*x + 3*c) + sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))
)*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 36*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(
3/2*d*x + 3/2*c)))^2 + 32*sin(3*d*x + 3*c)*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 16*s
in(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 8*cos(3*d*x + 3*c) + 1)*log(cos(1/3*arctan2(si
n(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2
- 2*sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1) - 48*cos(3/2*d*x + 3/2*c)*sin(3*d*x + 3*
c) + 80*cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))*sin(3*d*x + 3*c) + 48*cos(3*d*x + 3*c)*si
n(3/2*d*x + 3/2*c) - 4*(3*cos(3/2*d*x + 3/2*c) + 5*cos(7/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))
) - 3*cos(5/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 5*cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), c
os(3/2*d*x + 3/2*c))))*sin(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 20*(4*cos(3*d*x + 3*c) +
 6*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 4*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(
3/2*d*x + 3/2*c))) + 1)*sin(7/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 12*(4*cos(3*d*x + 3*c)
+ 6*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 4*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos
(3/2*d*x + 3/2*c))) + 1)*sin(5/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 24*(3*cos(3/2*d*x + 3/
2*c) - 5*cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), c
os(3/2*d*x + 3/2*c))) - 16*(3*cos(3/2*d*x + 3/2*c) - 5*cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2
*c))))*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 20*(4*cos(3*d*x + 3*c) + 1)*sin(1/3*arct
an2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 12*sin(3/2*d*x + 3/2*c))*B/((16*sqrt(2)*a^2*cos(3*d*x + 3*c
)^2 + sqrt(2)*a^2*cos(8/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 36*sqrt(2)*a^2*cos(4/3*arct
an2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 16*sqrt(2)*a^2*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(
3/2*d*x + 3/2*c)))^2 + 16*sqrt(2)*a^2*sin(3*d*x + 3*c)^2 + sqrt(2)*a^2*sin(8/3*arctan2(sin(3/2*d*x + 3/2*c), c
os(3/2*d*x + 3/2*c)))^2 + 36*sqrt(2)*a^2*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 32*s
qrt(2)*a^2*sin(3*d*x + 3*c)*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 16*sqrt(2)*a^2*sin(
2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 8*sqrt(2)*a^2*cos(3*d*x + 3*c) + sqrt(2)*a^2 + 2*
(4*sqrt(2)*a^2*cos(3*d*x + 3*c) + 6*sqrt(2)*a^2*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) +
 4*sqrt(2)*a^2*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + sqrt(2)*a^2)*cos(8/3*arctan2(sin
(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 12*(4*sqrt(2)*a^2*cos(3*d*x + 3*c) + 4*sqrt(2)*a^2*cos(2/3*arctan2
(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + sqrt(2)*a^2)*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x
 + 3/2*c))) + 8*(4*sqrt(2)*a^2*cos(3*d*x + 3*c) + sqrt(2)*a^2)*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d
*x + 3/2*c))) + 4*(2*sqrt(2)*a^2*sin(3*d*x + 3*c) + 3*sqrt(2)*a^2*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/
2*d*x + 3/2*c))) + 2*sqrt(2)*a^2*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(8/3*arctan2
(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 48*(sqrt(2)*a^2*sin(3*d*x + 3*c) + sqrt(2)*a^2*sin(2/3*arctan2
(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*s
qrt(a)))/d

Giac [F]

\[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\int { \frac {{\left (B \sec \left (d x + c\right ) + A\right )} \sqrt {\sec \left (d x + c\right )}}{{\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {5}{2}}} \,d x } \]

[In]

integrate((A+B*sec(d*x+c))*sec(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(5/2),x, algorithm="giac")

[Out]

integrate((B*sec(d*x + c) + A)*sqrt(sec(d*x + c))/(a*sec(d*x + c) + a)^(5/2), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {\sec (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{5/2}} \, dx=\int \frac {\left (A+\frac {B}{\cos \left (c+d\,x\right )}\right )\,\sqrt {\frac {1}{\cos \left (c+d\,x\right )}}}{{\left (a+\frac {a}{\cos \left (c+d\,x\right )}\right )}^{5/2}} \,d x \]

[In]

int(((A + B/cos(c + d*x))*(1/cos(c + d*x))^(1/2))/(a + a/cos(c + d*x))^(5/2),x)

[Out]

int(((A + B/cos(c + d*x))*(1/cos(c + d*x))^(1/2))/(a + a/cos(c + d*x))^(5/2), x)